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Abstract 

Tracy- Widom (TW) equations for one-matrix unitary ensembles (UE) (equivalent to 
a particular case of Schlesinger equations for isomonodromic deformations) are rewritten 
in a general form which allows one to derive all the lowest order equations (PDE) 
for spectral gap probabilities of UE without intermediate higher-order PDE. This is 
demonstrated on the example of Gaussian ensemble (GUE) for which all the third order 
PDE for gap probabilities are obtained explicitly. Moreover, there is a second order 
PDE for GUE probabilities in the case of more than one spectral endpoint. 

This approach allows to derive all PDE at once where possible, while in the method 
based on Hirota bilinear identities and Virasoro constraints starting with different bi- 
linear identities leads to different subsets of the full set of equations. 
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1 Introduction: Schlesinger-TW Equations and r-functions 

We consider the unitary ensembles (UE) of Hermitian random matrices. Various large n 
limits for this case belong to important widespread universality classes: bulk limit leads to 
equations for the sine kernel [El [15], [9] whereas edge limit (largest eigenvalues) leads to Airy 
kernel [HI US, [9]. One studies the operator K J with the kernel K J (x,y) = K(x,y)xj c (y), 
where K(x,y) = (ip(x)ijj(y) — ip(y)<p(x))/(x — y) and XJ c (y) is equal to zero on J and one on 
its complement J c in K. Important for this theory is the resolvent kernel of K (we will use 
short-hand notation K for K J ), R(x,y), the kernel of K{1 — K)' 1 , which is defined by the 
operator identity 

(I + R)(I-K)=I. (1) 
One introduces [EJ IS] the auxiliary functions 

Q{x- J) = (/ - X)-V(x), P(x; J) — (I — K)- 1 ^), (2) 
and auxiliary inner products [TH], which are functions of only the endpoints of J, 

U = (q, (j>X.Jc), v = (q, tpx.Jc) = (P, <J>X.Jc), w = (p, ^XJ-)- (3) 
Then the resolvent kernel R(x, y) is [8j [7J [15] 



, Q(x; J)P(y; J) - P(x; J)Q(y; J) . c „ , , m 

R[x, y) = , (x,y e J , x t y) (4) 

x-y 

and (prime denotes the derivative w.r.t. x, e.g. Q'(x; J) = dQ(x; J)/dx) 

R(x, x) = P(x; J)Q\x; J) - Q(x; J)P'(x; J). (5) 

The Tracy- Widom (TW) equations for one-matrix case are the equations for the functions of 
the endpoints — R(dj, a^), 

Qi = Q( a f, J )i Pj = p ( a f, J )i ( 6 ) 

and the auxiliary functions u, v, w. Among them there are universal equations valid for any 
Hermitian one-matrix unitarily invariant model. They read [15] 

dq 3 



da 



k 



da k 



-l) h R{aj,a k )q k , (7) 
-l) k R(aj,a k )p k , (8) 
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for j 7^ k and 



R jk = R(a j ,a k ) = qiPk m \ (9) 



-ly-'d, \nr J = R sj = R( aj ,a s ) = Pi ^-_ - q^, (10) 



3a, - <"W <»> 



<9w 

£ - (12) 

By definitions of the matrix integral over J n (or J°°), r J , and the resolvent operator R, we 
have the fundamental relation between them, expressed by the first equation in formula flTUJ). 
This is just a consequence of the defining connection, the equality of two different expressions 
for the probability of all eigenvalues to lie in J, 

r J 

det(I-K J ) = —, (14) 

T 

where r is the corresponding matrix integral over whole M n . The integrals r J and r can be 
considered as particular values of certain r-functions of integrable hierarchies, see e.g. [Tl|3l|T6]. 
The other equations are not universal, their particular form depends on the potential V(x). 

Let 

q' 3 = Q\a f J), T f j = P'{a j ;J), (15) 
Then the non-universal equations are 



dqj_ 
dctj 

dpj_ 
dcij 



q' j -J2(- 1 ) kR ( a ^ a ^, (16) 

k+j 

p;-^(-l) fe ^-,afeK, (17) 

k+j 



along with 



R(aj, dj) = pjq'j - qjp'j + ^(-l) /c i?(a j , a k )(qjp k - Pjq k ). 
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The above equations are equivalent to the Schlesinger equations involving 2x2 matrices Aj, 



see [TO] and [12] for details. 

The paper is organized as follows. In section two we show our general form of equations 
derived from the above. In section 3 we specify it to the GUE ensemble and obtain explicitly 
all independent third order PDE for its gap probabilities, for arbitrary number of spectral 
endpoints. In section 4 we consider the simplest representative example of two endpoints for 
GUE in more detail, show that some of the previous third order equations are redundant 
in this case and expose also a second order PDE (lower order than Painleve IV analogs!) 
satisfied by GUE probabilities. In section 5 we derive the results of section 2 which are the 
basis for all the rest. In the appendix we demonstrate the application of Adler-Shiota-van 
Moerbeke approach PQ [3] to GUE probabilities based on our system from [11]. The last 
sections presents the conclusions. 

2 The results. General case 

The system of Schlesinger-TW equations can be reduced in general to the following convenient 
form, which, at least for the Gaussian case, will immediately give all the 3-rd order PDE in the 
logarithm of the gap probability, In P, without intermediate higher-order PDE to integrate. 

Theorem 1. There are several series of universal equations independent of potential (i.e. of 
Xj, Yj defined below): the two u l-point" equations, 




(19) 



djudfr = GjdjjV - <)/!<), F. 



(20) 



G J 2 = (9 J F) 2 + 4F(^) 2 , 



(21) 



and the four series of "2-point" equations, 



P jt = djFdtF - GjGi + 4Fdjvd lV = 2F(a j - a^d^T, 



(22) 



GidjF - GjdiF = 2F( aj - a,)<9>, 



(23) 



djGi - diGj = 2(a,j - a^d^v, 



(24) 



(aj - ai)(djGi + diGj) = A{divd j F - djvdiF). 



(25) 
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Also there are two series of non-universal equations where the data dependent on the potential 
("spectral curve") enter, let Xj = q'j/qj, Yj = p'j/pj, then 

d J Vv = d 3 v(X J +Y J ), (26) 



(djBo - ajdjV) T = djviYj - Xj), 



(27) 



Here Bk = ^2jO>j dj, D — B-± — ^2jdj, T = lnr J , definitions of F and Gj are somewhat 
different for finite size n and infinite matrices. For finite n 



F = UW 



T n+l _ T n-1 
1~n 



Gj = WdjU - UdjW, 



u = u 



— (j i//\-l,„l /Tn+l r n-l f, T n+l/ T n+l\ 

Tn+lTn-1 ( T n-\l T n~\ 



w = w n = ^,{I-K J n )-^) 



1 , 



v = v n = (tp, (I - KiY^) = (/ - KiY V) = ~ In ^ 



t=0 



For infinite ensembles 



(28) 
(29) 

(30) 

(31) 

(32) 



F = —uw, 



(33) 



Gj = udjW — wdjU. (34) 

We observe that we have 4 series of 1-point equations (depending on one subscript j) and 4 
series of 2-point equations (depending on two subscripts, j and I). This means that for the case 
of iV spectral endpoints the total number of equations here is 4iV + 4iV(iV— 1)/2 = 2iV(iV + l). 
The model-dependent functions are supposed to be ultimately expressed in terms of the 
variables T, v, F and {Gj,j = 1, . . . ,N}. Then the number of dependent variables in the 
equations is going to be equal to iV + 3. We get only an upper bound on the number of 
independent PDE for the single function T only, assuming that the other N + 2 variables are 
eliminated, reducing the number of equations by N + 2: 

# equations for T < 2iV(iV + 1) - (N + 2) = 2N 2 + N - 2. 
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As we will see below on the example of GUE, the actual number of independent equations 
for T is less than this bound. It seems difficult to determine the exact number in general. 

Multiplying the equations ( 126]) and ff27|) by aj with an integer k > and summing over 
the endpoints aj, one gets 

Theorem 2. The model- specific equations for one-matrix UE can be presented in the following 
general form in terms of operators {B k , k > —1} of Virasoro subalgebra: 

B k ^Vv = $ fe , (35) 

(Bk-tBo - B k V)T = T k , (36) 
where we denoted $> k = ^ ■ a k jdjv{Xj + Yj) and Y k = J2j a> j9j v (Xj ~ Xj). 

Multiplying the universal "2-point" equations by a^a™, then summing over j and I and intro- 
ducing quantities 

Gk = 2j a j G i 

3 

leads to 

Theorem 3. The four common series of equations for all one-matrix UE in terms of operators 
{B k ,k > —1} of Virasoro subalgebra read: 

B k -iFB m -iF - G k G m + 4FB k -ivB m -iv = 2F{B k+1 B m ^ + B k ^B m+1 - 2B k B m )T, (37) 

G m B k -\F — G k B m _\F = 2F{B k B m ^\ — B k _\B m + £> fc+m _i)t>, (38) 
B k ^G m - B m -xG h + (k- m)G fc+m _i = 2(B k B m ^ 1 - B k ^B m + B fc+m _i)u, (39) 

Bm-\G k+ i—B m G k —G k+m = 2(B m -ivB k -iF—B k -ivB m -iF) — (B k+ iB m -i+B k -iB m+ i—2B k B m )v. 

(40) 

All the other equations that UE satisfy in terms of the operators B k only, are combinations of 
these equations and their derivatives. 
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It is worth to write out explicitly the components k, and k, k of the "multivariate 
Painleve-type" series (137)) . which contains indeed all multivariable generalizations of Painleve 
equations: 

Bk-iFVF - G k G + AFB k ^vVv = 2F(B k+1 V + B k ^B x - 2B k B )T, (41) 

(B k ^F) 2 -61 + AF(B k ^v) 2 = 2F(B k+1 B k ^ + B k ^B k+1 - 2{B k ) 2 )T. (42) 

One can eliminate Gj- variables using equations (12~TT) . (|22p and f[2"3"]) and obtain 2 indepen- 
dent series of universal PDE containing only T and functions v and F, which, for finite size 
ensembles, are related to T as [12] t> = — d t T\ t=0 and F = d 2 t T\ t=0 , where t is the first time of 
one- dimensional Toda lattice hierarchy [3] [11]. They can be written as e.g. 

Theorem 4. For any two spectral endpoints aj and ai, the following two PDE hold for all 
one-matrix UE: 

(d 2 t v) 2 - (a, - a{) 2 {dlTf + Ad^Tdjvdtv = 0, (43) 

(a, - a0 2 (F(d» 2 - d^TdjFdtF) - (divdjF - d.v^F) 2 = 0. (44) 
All the other such equations that UE satisfy are combinations of these and their derivatives. 

Remark. In the previous paper we used a higher-order PDE3, 

djV^F + 2(d J v) 2 ) = d 3 Fd 2 3 v - djTGj, (45) 

which is in fact the eq. (120]) where djGj has been eliminated using the derivative of eq. (|2T|) 
w.r.t. aj. Then from eqs. (|2"T|) and (145 p one obtains the higher-order universal "1-point" PDE 
of our previous paper [TJ] , 

(dAd^F + 2{d jV ) 2 ) - djFdjjv) 2 = (d,T) 2 ((d,F) 2 + AF{d 3 v)% (46) 

while the eqs. (143]) and (144]) are "2-point" equations. 

know it is higher-order due to the connection with Toda lattice since in general F = df t T\ t= Q, where 
t is the first Toda time. 
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3 Results for the Gaussian case 

For the Gaussian case the expressions for the non-universal quantities Xj and Yj can be 
found [15] and written as 

d jV (Xj + Yj) = 2djF, 

djv(Y j -X j ) = 2(G j + a j djv), 
so the two non-universal equations (1261) and ff27|) give, respectively, 

Vv = 2(F-n/2), 

2Gj = {djB - ajdjV) T - 2a j d j v. 
Then the linear equation (|24|) together with the last two relations result in formulas 

2v = VT, (47) 
4F = V 2 T + 2n, (48) 

2Gj = (djB - 2ajdjV) T. (49) 

Thus, all auxiliary variables are expressed in terms of T. Substituting these expressions into 
the other universal equations one finds the following two theorems: 

Theorem 5. The logarithm T of a spectral gap probability for GUE satisfies the following 
"1-point" 3rd order equations: 

djBoTd^VT - djVTd]jB Q T = djTdjV 2 T - 2(d j VT) 2 . (50) 

(djV 2 T) 2 - 4((<9 i £ - 2 aj djV) T) 2 + 4(V 2 T + 2n)(d J VT) 2 = 0. (51) 
All the other such equations it satisfies are combinations of these and their derivatives. 

Eq. fl50|) . however, becomes trivial in the case of one endpoint, while in the case of two 
endpoints the two such equations are a consequence of the other ones, see the next section. 
Whether this is the case in general, remains unclear at the moment. 
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Theorem 6. The logarithm T of a spectral gap probability for GUE satisfies the following 
"2-point" 3rd order equations: 

djV 2 Td{D 2 T - 4 (djB - 2a j d j V) T (diB - 2aAV) T + A(V 2 T + 2n)d j VTd{DT = 

= 8(Oj - a x ) 2 {V 2 T + 2n)d){T, (52) 

(diB - 2aAV) TdjV 2 T - (djB - 2a 5 dp) Td{D 2 T = 2(a i - a{){p 2 T + 2n)d 2 {DT, (53) 

2(a i - o/)(Sf,B - («j + ai)d 2 {D)T = d{DTdp 2 T - dpTd{D 2 T. (54) 
t4// i/ie other such equations it satisfies are combinations of these and their derivatives. 

We can now count the total number of third order PDE for the logarithm of GUE gap 
probabilities. We got two 1-point series and three 2-point series of equations, therefore in the 
case of N endpoints a, one has 2iV + 3iV(iV - l)/2 = iV(3iV + l)/2 independent PDE. This 
should be reduced by 1 for N = 1 since then (1501) is trivial and we are left with (15T|) only. 

The last three series of PDE have their counterparts in terms of i3-operators only. Since 
for the case of GUE we have 

$ fc = 2B k ^F = B k -iDv = B k _ 1 V 2 T/2, (55) 

T fc = 2(G k + B k v) = 2G k + B k VT = (B k -iB - B k V)T, (56) 

2G k = (B k ^B - 2B k V)T, (57) 
we obtain immediately, either from theorem 3 or from theorem 6, 

Theorem 7. The logarithm T of a spectral gap probability for GUE satisfies the following 
equations in terms of operators {B k ,k > —1} of Virasoro subalgebra: 

B k ^ 2 TB m ^ 2 T-A(B k ^B -2B k V)T(B m ^B -2B m V)T+A(V 2 T+2n)B k _PTB 

= 8(V 2 T + 2n)(B k+1 B m ^ + B k ^B m+1 - 2B k B m )T, (58) 
(B m ^B - 2B m V)TB k _ 1 V 2 T - (B k ^B - 2B k V)TB m _ 1 V 2 T = 
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= 2(V 2 T + 2n)(B k B m ^ - B k ^B m + B k+m ^)VT, (59) 
(B m ^(B k B - 2B k+l V) - B m (B k -iB - 2B k V) + B k (B m ^B - 2B m V)- 

-B k ^(B m B - 2B m+1 V)) T = B m ^VTB k ^V 2 T - B k ^VTB m ^V 2 T, (60) 

All the other equations in terms of the operators B k it satisfies are combinations of these and 
their derivatives. 

In fact, eq. (12 ip can be considered as eq. ( I22p for CLj — Qj\ , therefore, summing ( 12 2 j) over 
both subscripts j and /, we get the component in a sense corresponding to (12"TI) also, namely, 
for GUE - the "diagonal" equation (|58|) with m = k. The total number of algebraically 
independent equations in the three series is thus 3N(N — l)/2 + N = N(3N — l)/2 in the 
case of N endpoints. It is less than the total number of 1-point and 2-point equations above 
by N, as should be since equations floTI]) do not have a counterpart here. 
E.g. the simplest members of the three series are, respectively, 

(£> 3 T) 2 - 4(£>T - B VT) 2 + A(V 2 T + 2n)(VT) 2 = 1Q(V 2 T + 2n)(B 1 V - (B ) 2 + B )T, (61) 



(VT - B VT)B V 2 T - {{B ) 2 T - 2B 1 VT)V 3 T = 2{V 2 T + 2n)(B 1 V - (B ) 2 + B )VT, (62) 



(ViBiBo - 2B 2 V) - 2-B ((-B ) 2 - 2B{D) + B\{V - B V)) T = V 2 TB V 2 T - B VTV 3 T. 

(63) 

Equation (|6T|) was first found by comparison of two 4th order PDE - the "boundary-KP" 
equation of [TJ and another one derived by the author from "boundary- Toda-AKNS" system 
of [UJ, see Appendix. Now the general procedure to obtain all the lowest order equations at 
once, wherever possible, is outlined. 

4 Example: Gaussian case with two endpoints 

We start with general equations of theorem 1. However, let us introduce notations which will 
be convenient for Gaussian two-endpoint case. If the two endpoints are aj and ai, then let 

£ + = aj + a h = aj -a h V = dj + d h £>_ = d s - d h 
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G-\- — Gj -\~ Gi, G— — Gj — G{. 

Then we have 

dj = 1/2(1? + ZL), di = \/2{V-VJ), Gj = 1/2(G+ + G_), G, = 1/2(G+ - G_). 



In the new notation equations of theorem 1 read as follows: 

(G + + G_) 2 = (£>F + XLF) 2 + AF(Vv + P_<u) 2 , $HD 

(G + + G_) 2 = (£>F + XLF) 2 + 4F(£>u + V^vf, (64) 

(£>F) 2 - (ZLF) 2 -Gl + G 2 _+ AF((Vvf - (£>-^) 2 )) = 2F^(V 2 T - X? 2 T), (|22l) 

G+V_F- GJDF = £-F{V 2 v -V 2 _v) } (J2D 

T?_G+-£>G_ = £_(X?\-£> 2 t>), (TjSD 

£-(£>G+ - T?_G_) = 4(VvV_F - V^vVF), (TJSD 



(2X; + £Lt;)(X? + £>_)(G + + G_) = (G+ + G_)(£> + £>_) 2 w -2(VT + V_T)(VF + V_F), (TJUD 

(£>u -£>_?;) (£>-£>_)(G+-G_) = (G + -G_){V -V_) 2 v - 2{VT -V_T){VF -V_F). (65) 
Adding and subtracting the first two equations above, one gets, respectively, 



G 2 + + G 2 _ = {VFf + (£>„F) 2 + AF((Vv + (V_vf), (66) 

G + G- = VFV^F + AFVvV^v. (67) 
Adding and subtracting (1221) from fl66l) gives, respectively, 

G+ = {VFf + 4F(£>i;) 2 - ^F(P 2 - V 2 _ )T, (68) 
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G 2 _ = (£LF) 2 + AF{V„v) 2 + ^F(V 2 - V 2 _)T. (69) 
The non-universal equations for the Gaussian two-endpoint case give 

4G+ = 2VT - £ + V 2 T - i_VV_T, (70) 

AG- = 2T?_T - i_V 2 T - i+VV_T - (V 2 T - V 2 _T) , (71) 

and already known expressions for v and F in terms of T. One sees again that (1241) is satisfied. 
Denoting also 

r = VT, S = VV_T, r_ = T?_T, A = V 2 T - V 2 _T, 
one finds expressions 

v = r/2, F = l/4(Vr + 2n), 
Vr - V_r = A, X?r_ = X>_r = S, 
Vv = Vr/2, V^v = S/2, (V 2 - V 2 _)v = 1 /2VA, 
V_S = V 2 r -VA, V 2 _v = l/2(Vr-VA), 
VF = l/AV 2 r, XLF = l/AVS, 
4G+ = 2r - i + Vr - ^S, 4G_ = 2r_ - £_T?r - £ + S - ^A, 
4VG + = -i+V 2 r - £_VS, 4VG_ = -£_V 2 r - £ + VS - i-VA, 

4V_G + = -t+VS - £_(P 2 r - VA), 4ZLGL = -£+(£> 2 r - VA) - i_VS - i_V_A - 4A 
Using all these expressions, we can rewrite our system of PDE as 

(2r - i+Vr - i_S) 2 = (V 2 r) 2 + A{Vr) 2 (Vr + 2n) - 4£!(£>r + 2n)A, (EHD 

(2r_ - i_Vr - i+S - ^A) 2 = (VS) 2 + AS 2 (Vr + 2n) + 4£ 2 (Vr + 2n)A, (PJ) 
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(2r - i + Vr - £-S)(2r_ - £_T?r - £ + S - ^A) = V 2 rVS + AVrS(Vr + 2n), (EZD 

(2r - i + Vr - i_S)VS - (2r_ - £_X?r - - i_A)V 2 r = 2£_(£>r + 2n)VA 1 (1231 

2(S£> 2 r - VrVS) = £-(£+VA - £-V_A - AA), (J25D 

and equations fl20|) and fl65l) reduce after a long tedious calculation to the following two 
equations, 

2^AV 2 r = {i-{Vr + A) - 2r_)£>A + S(£_£>_A + 4A), (72) 
2f_AX>S = (^S - 2r)VA + Z>r(f_£LA + AA). (73) 

Let us denote also 

F = Vr + 2n, D_ = £_X?_A + 4A 
Then, expressing PA from (1231) 

2i_FVA = G+VS - G_V\ ([23]) 

and then D_ - from ( 1231) . 

2^FD_ = (£ + G + + AFVr)VS - (£+G_ + 4FS)X? 2 r, (1231) 

Substituting these expressions into (172|) and (173|) . one sees that the last equations become 
trivial. Thus, they are also redundant, i.e. are consequences of the rest of the system. 

Remark. One can notice a lot of similarity with the final system of PDE for joint largest 
eigenvalue probabilities of two coupled GUE, see the Corollary in section 6 of [T3] ■ This is to 
be expected in light of the representation of the matrix model with one finite size spectral 
gap in terms of models of two coupled matrices in [5], p. 6745. More, however, remains to be 
understood about these relations. 

We are left with five independent PDE involving four independent third derivatives of T. 
This means that there should be one second order PDE in T, i.e. equation of lower order 
than multidimensional analogs of Painleve IV. Indeed, this equation can be easily obtained, 
if we express (V 2 rVS) 2 from eq. (I67p and compare with its expression obtained from the 
corresponding product of (16"H1) and (15^1) . It reads as 
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(VrG. - SG+) 2 = £A((GL) 2 - (G+) 2 + AF((Vrf - S 2 - (74) 

It becomes trivial upon the identification of the two endpoints as it should be. The analogous 
equation for GUE can be obtained also in general, for any number of endpoints. This follows 
by expressing (djFdiF) 2 from ( )22|) and comparing with its expression from the corresponding 
product of (12TI) and (164")) . now for any pair of endpoints cij and a;. Thus, we obtain 

Theorem 8. The logarithm of a GUE gap probability T satisfies a second order PDE, 
{d{DT{djB T - lafipT) - d/DT{dfi Q T - 2 ai d{DT)) 2 = 



4(a-a z ) 2 <9j z T(-(c> i £oT-2a^ 

(75) 

5 Transformation of the TW-Schlesinger system 

Here mostly the derivations of the statements in the previous sections are gathered. 
An important identity is obvious from definitions in section 1: 

(djv) 2 = djudjW. 

Taking the second derivatives of u, v and w and using the equations for derivatives of qj and 
Pj, one can obtain: 

s^l n 1 s~J n « s~l n t s~i n « 



2 J 

aj - 


- ai 


djudiw — 


diudjiu 


a j ~ 


a x 


^djvdiw - 


- divdjW 


dj - 


- ai 


(-iy<&5j 









d 2 3 w = 2(-iy Pj p' J -J2d>- 
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The fisrt derivatives of T can be written in two different ways: on the one hand, by fllOl) . (1161) 
and (PU), 

d T = ( \y II,, = (-ly-'ip^ - qfltj) - £(-l) i+, i§(a,- - o,), 
on the other hand, by ([TDD and AHMED, 



djT = ( \y '% = ( I)' (p/h'lj - qjdM) 



djudjjW - djwdjjU 84 
2d~v ~ 



2 3 



d jU 



The terms in the sum of the first formula can be also represented as 

djudiiu + diudjW — 2djvdiv 



[-iy + %( aj - ai ) 



a,-, — a\ 



Switching to the variables U, W, Xj, Yj introduced above, one can rewrite the equations as 
follows: 

djudiw = -djUdiW, (d jV ) 2 = -djUdjW, 



n d 3 Udiv 


— diUdjV 


aj 


- ai 


djWdtU - 


- diUdjW 




- ai 


2 3 


- divdjW 



a,j — ai 



d] j U = 2d j UX j -J29- i U, 



dj J v = d j v(X j + Y j )-Y,d 



8% j W = 2d j WY j -Y,dflW, 
djT = djviYj - X 3 ) - - a,) 

djW&JJ - diUdlW d 



Q T = ^3- -33- = ln 

3 2d,v 2 3 



d 5 U 



d 3 W 
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Now we express everything in terms of new functions F and G 3 defined in the statement of 
theorem 1. The first simple identity becomes 



(d jV 



which is nothing but the equation (T2~Tj) . Mixed second derivative equations now read as 



d 2 rr _ + Gj)div - (diF + Gj)djV 

31 (ai-ai)F 



2 _ GidjF - GjdiF 
31 2(a i -a l )F 



31 (o,--o,)F 



Since 



wd%u - udlw = + diGj), wd^u + u%w = d) x F - d3FdlF 2F ^ 

it follows that 

djGi + d lGj = A 8 - 3FdlV ~ 9lFd3V 



CLj — &\ 

, , . djFdiF - GjGi GAv - GidjV 

dlF - 3 — = 2^^ 76 

2F dj - a t 



By definitions, either (125]) and ( 129|) or ( 133|) and ( 134)) . there are also equations of Adler-van 
Moerbeke type [I], 

a, (%) = a, fay (77) 



see [13] for the discussion of their analog for coupled GUE joint distributions first derived 
in jl]. The last equation together with (123 j) means that a general linear PDE holds: 

9,-G, - SkGj = 2( Qj - a,)Sf,v. d21 
The diagonal equations then read: 

d 3 G 3 = (d jF + Gj)Xj - - G,)^- - 2 £ (78) 
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W - = W + ^ + (%F - G,« - 2 £ G '»'-^" , (79) 

%v = dAx 3 + is) - E 2{a .- ai)F ■ 

Using the above equation ( 1231) and the diagonal equation for djjV, the second derivatives d^v 
can be summed over all I, giving 

Vd j v = d j v(X j + Y j ). flUD 

For djT one has expressions 

m m 2{a > ~ ai)F 

(80) 

and 



2Fd j Fd j G j - G^Fd^F - («9,F) 2 + G}) d 3 Fd 3 G 3 - G^F + 2{d,vf) 

AFdjv 

If we sum over j the first of these two equations, the last sum on its right-hand side cancels 
out since its terms are antisymmetric in indices j and I, so we obtain 

VT = y £ i d j v(Y j -X j ). (81) 

3 

The second derivatives of F, entering eqs. (1761) and (1751) . can in fact be removed from con- 
sideration. This is because we can take the corresponding derivatives of equation (J2T1) and 
compare with (1761) and ( 1791) . The result is that equation ( 1761) is redundant - it follows from the 
derivative of (12ip w.r.t. a; and the other equations. As for the equation (1791) . its comparison 
with the derivative of (l2Tj) w.r.t. a,-, 

2d j Fd] J F - 2G j d J G j + AdjF(d jV ) 2 + SFdjV^v = 0, 

gives immediately the equation (120]) . which is thus also proved. As we will see, at least for 
Gaussian case using eq. (1201) is preferable since it encodes a third order PDE for T, while (1791) 
corresponds to a fourth order PDE in T. So we keep (1211 and then (1791) is also redundant. 
Thus, what remains to prove theorem 1 is to derive equations (j2"2"|) and (J2"7|) . 



8F 2 6»t> 
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Proof. The simplest way to prove (1221) and (127)) is to differentiate (1101) over a/, I ^ j, and then 
use fl7J), (JS]), which gives 

djiT = (- ,<)],!> j - !>,<)],(/, + dtfjdjPj - dipjdjqj) = 

= (-^^(QjdjiRjiPi) - Pjdj(Rjiqi) + RjiiqidjPj - pidjqj)), 
but ([TJ and (jSJ) also imply that 

- Pjdjqi = 0, 

and therefore, differentiating Q, we get 

^ ^ dj — (Xi CLj — 0,1 

Applying the last two equations results in the fundamental relation 

d\T = -(-iy +l R% (82) 

which is equivalent to (122)) if we express Rji in terms of F, Gj, G[ and v. Now, substituting 
d 2 {F into the summand of (1801) and gathering terms with derivatives of T together, one obtains 
equation (127)) . □ 

It remains to prove theorem 4, since all the other results are consequences of theorem 1. 

Proof. We first multiply eq. (122)) by djFdiF, take square of eq. (123)) and then compare them, 
getting 



2G j d l F ■ GidjF = 2(d j Fd l F) 2 + 8FdjvdiF ■ d^djF - 4F( aj - aifd^TdjFdiF = 

= G 2 (d t F) 2 + Gf(djF) 2 - AF\a, - a,) 2 (9» 2 . 

Then we use eq. ( 12T)) to express G 2 and G 2 , after which some terms cancel out and the last 
equality above becomes eq. ( 144")) of theorem 4. 

To derive eq. (|4"3~)) . we first write down a direct consequence of (122)) . 

(djFdiF - GjGi) 2 = AF 2 ((a, - ai ) 2 d 2 {T - 2d J vd l v) 2 , 
then use the square of ( 125)) again, obtaining 
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Applying again eq. ( 121 j) to eliminate G| and G z 2 from the last equation leads to eq. (143 j) after 
cancellation of several terms. 

It is clear that the number of such independent equations is two since we used four equa- 
tions algebraically eliminating two variables - Gj and Gi from them. 

□ 

6 Conclusions 

All independent lowest order PDE for GUE spectral gap probabilities are explicitly obtained. 
This gives also an example of explicit integrable hierarchies in terms of isomonodromic "times" 
- the spectral endpoint variables. 
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7 Appendix: PDE for Gaussian matrices from bilinear 
identities and Virasoro constraints 

In [11] from Toda-AKNS system we obtained the following system of three equations for 
Gaussian one-matrix gap probabilities: 

B 2 _ 1 hxri = 2{2U n W n -n), (83) 
E\U n = (-2B + An)U n - W 2 n W n , (84) 

B\W n = (2B + 4n)W n - 8U n W 2 , (85) 

where U n = t^ +1 /t^, W n = t^_ x /t^, and the operators Bk contain partial derivatives w.r.t. 
the boundary points of the intervals: 

i(^{endpoints} 

It can be solved in general. Let F = UW = U n W n , G = WB^U - UB-xW, G = 
WBqU — UB W, T = lnr^. We will need the commutation relation: 

B B^ = B_x(B - 1). (86) 
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The combination W ■ flM|) — U • (EST) together with (155)) gives 

S-iG = -2B F = - l -B,B\T = -^B^(B - 1)B^T, 
or, taking into account vanishing of all the derivatives at infinity, 

G = ~(B - l)H_iT. (87) 
The combination W ■ ( 18~4l) + £/ • ( 1851) and the above definitions of new variables lead to 



B\F - ° 2 = -2G - 16F(F - £). (88) 



IF ' 2 

Another combination, B-{W ■ (JH3J) + B-iU • (JHS]), results in 



— ) = -2(B_ 1 WB C/ - B^UBoW) - 8(F - \)B- X F. (89) 
If we try to use the expression 

on the right-hand side of (1891) . we get just eq. (IHHj) back again and nothing new. One needs 
to use another expression for the right-hand side of (1891) : 

2(B- 1 WB U - B-xUBoW) = B^G - (B + 1)G, 

which can be obtained by moving the differentiation operators in B-iWBqU — B-\UBqW 
to the left in two possible ways and then matching them and using eq. ( 186|) . Then eq. (189]) 
becomes 

2£_iG = 2(B + l)G - 88_i(F - '-f - £_i 1 i - B - lF ^ ~ G * 



2 ' \ 2F 

Luckily, here (Bo + 1)G can also be expressed as £>_i-derivative, using eqs. ( 1H7|) and ( 1861) : 

2(£ + 1)G = -(£ 2 - l)B-{T = -B_!((B - I) 2 - 1)T = S-i(28o - 
so eq. ( 189|) is a total derivative and its integration gives 

2G = (2i3 - Bl)T - 8(F - \f - ' °' - . (90) 

Finally, substituting eq. ( l90l) into eq. ( 1881) and using eqs. ( 1831) and ( 1871) . we obtain a single 
equation for T = lnr^: 
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{B\ - ABl + 8B + MB\)T + 2{B\T) 2 - ( ^ lT) J^l^'^ = °- ( 91 ) 

This is a fourth order PDE as is the "boundary-KP" equation of Adler, Shiota and van 
Moerbeke [I], 

(B\ + SnB 2 ^ + 12B 2 + 2AB - 16£_i#i) lnr„ J + Q(B\Tf = 0, (92) 

derived from the KP equation. Comparing the eqs. (1911) and 0921) . we obtain a third order 
equation for T. Using another commutation relation: 

B-iBi = BiB-x + 2B , (93) 

it can be written as 

{B\Tf - 4((fib - l)B^Tf + A(B\T + 2n)((B 2 _ 1 T) 2 + (B 2 -B - B X B_ X )T) = 0. (94) 

This last equation is a direct multidimensional analog of the well-known Painleve IV equation 
for single endpoint Gaussian case [151 [1] : 

(r»y _ 4 (£ r ' _ r )2 + 4 ( r ')2( r ' + 2n ) = o, (95) 

where r = T' and 'prime' means the derivative w.r.t. the endpoint £, since for the case of 
single endpoint it turns into the equation (|95|) . The last term in the parentheses of eq. (|94|) . 
which has no analog in eq. (195]) . turns then to zero. 

However, the equation (|94|) is not the only third-order PDE satisfied by T, see the others 
for the 2-endpoint case in the sections above. Thus, this situation that the gap probabilities 
for more than one spectral endpoint satisfy several different PDE is very general and applies 
for both single and coupled ensembles. 
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